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Abstract. We address the problem of constructing numerical integrators for non- 
holonomic Lagrangian systems that enjoy appropriate discrete versions of the geomet- 
ric properties of the continuous flow, including the preservation of energy. Building 
on previous work on time-dependent discrete mechanics, our approach is based on a 
discrete version of the Lagrange-d'Alembert principle for nonautonomous systems. 

1. Introduction. In the last years Geometric Integration has grown to be a very 
large and active area of research, with a rich variety of approaches taken and top- 
ics covered |51 128) . Among the various viewpoints, the variational integrators ap- 
proach has revealed to be very powerful 21 . This point of view is not confined 
to Lagrangian and Hamiltonian (conservative) systems, but also admit extensions 
to multisymplectic geometry and PDEs, as well as to systems subject to external 
forces and dissipation (see [231 f° r a recent overview on the subject). 

The treatment of problems with constraints has also been an important issue in 
the area. Holonomic constraints have received a great deal of attention |1(J1 1121 1161 
IT7I l2l)l I5U| , motivated by their presence in applications such as molecular dynamics 
and planetary motions. The treatment of nonholonomic constraints has also been 
in the agenda of the Geometric Integration community (see, for instance, |241 133| ). 
Following the variational approach to discrete mechanics, we proposed in [H] a class 
of nonholonomic numerical integrators enjoying discrete versions of some of the 
geometric properties of the continuous flow. These include the evolution of the 
symplectic form along the flow, and the fulfillment of a discrete version of the non- 
holonomic momentum equation [3], which in the case of horizontal symmetries gives 
rise to conservation laws. However, these integrators do not preserve the energy, 
which is a natural conserved quantity of the continuous flow. This is not surpris- 
ing, since (fixed time-step) variational integrators themselves do not preserve the 
energy either. A different approach based on the technique of generating functions 
is proposed in [T5| . 

In this paper, we address the problem of energy conservation building on pre- 
vious derivations on time-dependent discrete mechanics and extended variational 
integrators [I'M 1181 12,'{j . Our main contribution is the construction of extended non- 
holonomic integrators derived from a discrete version of the Lagrange-d'Alembert 
principle for nonautonomous systems. We focus on investigating the relationships 
between the discrete and the continuous mechanics. The special feature of these 
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integrators is that, in addition to inheriting good properties with respect to the 
symplectic form and the nonholonomic momentum, they also preserve the energy. 

The paper is organized as follows. Section gives a brief introduction to me- 
chanical integrators for unconstrained systems and explains the necessity of allowing 
variable time steps to design algorithms which preserve at the same time the sym- 
plectic form, the momentum and the energy. The basic theory on time-dependent 
variational integrators is also presented. In Section|31 we propose a discrete version 
of the Lagrange-d'Alembert principle for nonautonomous constrained Lagrangian 
systems. This principle leads us naturally to the extended discrete Lagrange- 
d'Alembert equations, which we term nonholonomic integrators. Section 0] presents 
an account of the geometric properties of these integrators, paying special attention 
to the energy conservation. Finally, Section [5] gives some concluding remarks. 

2. Mechanical integrators. In this section, we briefly introduce some common 
notions and results from the literature on Geometric Integration. For further refer- 
ence, the reader is referred to [111 121)1 l2*2l I29| . Given a symplectic manifold (P, u>) 
and a Hamiltonian function H : P — ► R, an algorithm Fh : P — > P, h 6 [0, ho), 
is called a symplectic integrator if each Fh : P — > P is a symplectic map; an en- 
ergy integrator if H o Fh = H; and a momentum integrator if J o Fh = J, where 
J : P — > g* is the momentum map associated with the action of a Lie group G on 
P. An algorithm having any of these properties is called a mechanical integrator. 

The choice of a specific integrator depends on the concrete problem under con- 
sideration. For instance, in molecular dynamics simulation, the preservation of the 
symplectic form is important for long time runs, since otherwise one may obtain 
totally inconsistent solutions. On the other hand, the exact conservation of momen- 
tum first integrals is essential to problems in attitude control in satellite dynamics, 
since this is the basic physical principle driving the reorientation of the system. 
However, one is in general prevented from finding integrators which preserve the 
three elements at the same time due to the following result. 

Theorem 1 (H3). Consider a Hamiltonian system with a symmetry group G such 
that the dynamics Xh is nonintegrable on the reduced space (in the sense that any 
other conserved quantity is functionally dependent on H). Assume that a numer- 
ical integrator for this system is energy- symplectic-momentum preserving and G- 
equivariant. Then, the integrator gives the exact solution of the problem up to a 
time reparameterization. 

Roughly speaking, this result means that obtaining a fixed time step energy- 
symplectic-momentum integrator is the same as exactly obtaining the continuous 
flow. This theoretical obstruction can be overcome by allowing for varying time 
steps as we will review below. 

2.1. Variational integrators. Mechanical integrators derived from discrete me- 
chanics have their origin in the works by Lee, Veselov and others (see |15M25ll^Tll32| 

and references therein). In the last years, they have been intensively studied and 
further developed to deal with more general situations We briefly 

review here the main ideas of this approach. A complete exposition can be found 
in the recent overview . For the sake of conciseness, we directly go to the time- 
dependent case, without presenting the autonomous situation. 

Let Q be an n-dimensional manifold, and consider the extended configuration 
manifold Q — M x Q. The extended discrete Lagrangian state space is Q x Q, 
with canonical projections 7Tj : Q x Q — > Q, i = 1,2. An extended discrete path 
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is a sequence of points in Q, i.e. a map c : {0, . . . , N} — > Q. We denote c(fc) = 
{tkilk) & Q, k = Q,...,N. Given a discrete path, the associated discrete curve is 
q : {to, . . . , tj^} — > Q, g(£fe) = q k . The extended discrete path space is defined by 

C d = {c : {0, . . . , N} -> Q | t fc+ i > t fc , ft = 0, . . . , N - 1} . 

The tangent space T c Cd to Cd at c is the set of all maps 5c : {0, . . . , TV} — > TQ such 
that Tq o 5c = c, where t-q : TQ — > Q denotes the canonical projection. Consider 
the space (Q x Q) 2 =QxQxQxQ with projections at : (Q x Q) 2 — > Q x Q, 
i = 1,2. The extended discrete second-order manifold of (Q x Q) 2 is defined by 
Qd = {w e (Q x Q) 2 | 7T2 o ai(w) = 7Ti o er 2 (w)}. Otherwise said, Qd is the set of 
points w in (Q x Q) 2 of the form w = (to, qo, h,qi, t\, qi,t%, 92)- 

An extended discrete Lagrangian system is given by a map Ld : Q x Q — > R. The 
extended action sum S : Cd — > R is then defined by, 
w-i JV-l 
S(c) = L d (c(k),c(k + 1)) = X! L d (t k ,q k ,t k+1 ,q k+1 ) . (1) 

fc=0 fe=0 

Theorem 2 ( 23 ). Given a C k extended discrete Lagrangian Ld '■ Q x Q — > R, 
> 1, t/iere exists a unique C k ~ 1 mapping DuEhLd ■ Qd — > T*Q and unique C k ~ 1 
one-forms < dj ld and Oj d on QxQ such that, for all variations 8c £ T c Cd of c € Cd, 

AT-l 

d«S(c) '(5c= ADELAi(4-i,9fc-i,*fe,9ft,*ft+i,9fc+i) ■ (St k ,Sq k ) 

k=l 

+ Qi~ d (tN-i,qN-i,tN,qN) ■ (8tN-i,5qN-i,5t N ,Sq N ) 

- Z d (*o, 9o,*i, 9i ) ■ (£to,£go,<fti,<fyi) ■ (2) 

The map -DDEL^d is called the extended discrete Euler- Lagrange map and the 
one-forms 0j d and ©£ d are the extended discrete Lagrangian one-forms. Locally, 

DDEhLd(t k -i,q k -i,t k ,q k ,t k+1 ,q k+ i) = [D 4 Ld(t k -i , qk-1 , tk, Qk)+ 

D2Ld(t k ,q k ,t k+ i,q k+ i)]dq k + [D 3 Ld(t k ^i,q k -i,t k ,q k )+DiL d (t k ,q k ,t k+ i,q k+ i)]dt k , 

Qt d (tk,qk,t k +i,q k +i) = D4Ld(t k ,q k , t k+ i,q k+1 )dq k+ i+D 3 Ld(t k ,q k ,t k+ i,q k+ i)dt k+ i, 

®~i d (t kl q k ,t k+ i,q k+ i) = —D2Ld(t, k , q k ,t k+ i, q k+ i)dq k — D\Ld(t k , q k , tk+i, q k +i)dt k , 

where Di denotes the differential with respect to the ?th variable, i — 1, . . . , 4. 

To ease the exposition, along the paper we consider smooth discrete Lagrangians. 
Discrete Hamilton principle. The discrete variational principle states that, 
given fixed end points (io,9o), (^iv,9Jv), the evolution equations extremize S. 

Otherwise said, we seek discrete paths c £ Cd which are critical points of the 
discrete action, dS(c) ■ Sc = for all variations 8c S T c Cd with 5c(0) = = Sc(N). 
From Theorem |21 we get the extended discrete Euler-Lagrange (EDEL) equations, 

DxmhLd(tk-i,qk-i,tk,qk,tk+i,qk+i) = 0, 1 < k < N - 1 , (3) 

which can be equivalently written us 

D 2 L d (t kl q k ,t k+1 ,q k+ i) + D 4 L d (t k -i, q k -i, t k , q k ) = 0, (4) 

DiL d {t k ,q k ,t k+ i,q k+ i) + D 3 L d (t k -i, q k -l, t k , q k ) = . (5) 

If we define the discrete energies of the system to be 

Et d (tk,qk,t k+ i,q k+1 ) = -D 3 Ld(t k ,q k ,t k+1 ,q k+ i) , 

^L d (^k,qk,t k+ i,q k+ i) = DiLd(t k ,q k ,t k+ i,q k+ i) , 
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then equation (JSJ) can be simply written as 

E Ld {tk-itqk-i>tkt<lk) = E£ d (t k , q k , t k+ i, q k+1 ) , 

which reflects the evolution of the discrete energies. 

Under appropriate regularity conditions on the discrete Lagrangian (see 23 ), 
the DEL equations induce an extended discrete Lagrangian map $>:QxQ^QxQ, 
(tk-i,qk-i,tk,Qk) >-> (tk,qk, t k +i,q k +i). The basic geometric properties concerning 
extended variational integrators derived from the EDEL equations are the following, 

Symplecticity: consider the restricted discrete action S : Q x Q — > R, 

S(t ,q ,h,qi) = S(c) , 

where c G Cd is the unique solution of the EDEL equations satisfying c(0) = (to, qo), 
c(l) = (tx,qx). From Theorem^] we compute dS = ($ N ~ 1 )*Oj jd - 6£ , and then 

where ilL d is the extended discrete Lagrangian two- form, Qi d = — ^©L d = — ^®L d - 
Therefore, extended variational integrators are symplectic |13U18ll25| . 

Extended Noether's theorem: Let G be a Lie group acting on Q, i/j : G x Q — > 
Q, and consider its diagonal extension to Q x Q, 

# : GxQxQ — >QxQ 

(g,t Q ,q ,t 1 ,q 1 ) i — ► (ip(g,t Q ,q ),tp(g,t 1 ,q 1 )) . 

The discrete Lagrangian is G -invariant if Ld(^f(g, to, qo, ti, qi)) — Ld{to, qo,ti,qx), 
for all g £ G, (t , q ), (ti, gi) € Q. The discrete Lagrangian is infinitesimally in- 
variant if %xq) =0 >^GB,where^ x g(to,g'o,*ij9i) = (£q(*o,«)),£q(*i,?i)) 
denotes the infinitesimal generator of \& associated with £. Clearly, invariant La- 
grangians are also infinitesimally invariant. Using dLd = Oj d — , one sees that 
an infinitesimally invariant Lagrangian defines a canonical discrete momentum map, 

Jh d - QxQ — > g* 

(to,qo,ti,qi) ' — > JL d (to,qo,ti,qi) : -> K 

£ ^ e L d ' %xQ = e L d ' %xQ ■ 

If is G-invariant, then it can be easily seen that \t*G]? = . This implies 
that 3d is Ad-equi variant. A second fundamental fact is that extended variational 
integrators preserve momentum [131 1181 l2~3"] . i.e. Jl a o $ = J^ d . 

Energy conservation for autonomous discrete Lagrangians: a discrete 
Lagrangian is called autonomous if it is invariant with respect to the additive 
action of R on the time component of Q, i/j : R x Q — > Q, ip(s,(t,q)) = (s + 
t,q). The associated discrete momentum map is given by Ji d (to, qo, t\, qi) — 
— E~£ (to, qo, ti, q\)dt\ = —E^(to,qo,ti,q±)dto- Noether's theorem thus gives 

^L d (^k,qk,tk+l,qk+l) = E Ld (tk-l,qk-l,tk,qk) ! 

or equivalently, E^ d {t k , q k ,t k+ll q k+1 ) = E^ d (t k -i,qk-i,tk,qk), i-e. the discrete 
energy is conserved by the extended variational integrators derived from an au- 
tonomous Lagrangian |13l 1181 12*3*] . 

3. A discrete Lagrange-d'Alembert principle for nonautonomous sys- 
tems. In this section, we propose a discrete version of the Lagrange-d'Alembert 
principle for nonautonomous discrete systems. We start by defining what we un- 
derstand by an extended discrete nonholonomic Lagrangian system, 
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Definition 1. An extended discrete nonholonomic Lagrangian system is a triple 
{Ld,T>d,T>), where Ld : Q x Q — > R is the discrete Lagrangian, I'd C Q x Q is the 
discrete constraint space and T> is the constraint distribution on Q. In addition, T>d 
has the same dimension as V and is such that (t, q, t, q) G T>d for all (£, q) G Q. 

Notice that the unconstrained discrete mechanics (cf. Section l^"T|l can also be 
seen within this framework, where V = TQ and T>d = Q x Q. 

Remark 1. The motivation for this notion of extended discrete nonholonomic La- 
grangian system is the following. When dealing with unconstrained systems, given 
fixed end points (to, qo), (ijv, Qn)j one extremizes the action sum S with respect to 
all possible discrete paths. This means that at each point (t, q) £ Q, the allowed 
variations are the whole tangent space T( t , q )Q- However, in the nonholonomic case, 
one must restrict the allowed variations at each point: these will be exactly given by 
the distribution of feasible velocities T>. On the other hand, the discrete constraint 
space T>d will impose certain constraints on the solution sequence {(tk,Qk)}- 

Here, we will only consider constraints which do not impose conditions on the 
time velocities, i.e. t*(2?) = TR, where T:Q = RxQ^Ris the projection onto 
the first factor, although most of the discussion can be also carried out in broader 
terms. The constrained discrete path space is the set of extended discrete paths 
which verify the discrete constraints, 

C d = {c e C d | c(k) eV d , 0<k<N}, 

and the set of allowed variations is given by 

V d = {5c e TC d | Sc(k) G V c{k) , < k < N} . 

Discrete Lagrange-d'Alembert principle. Given fixed end points (to,9o) and 
(ijV)?Jv) ; the discrete Lagrange-d'Alembert principle consists of cxtrcmizing the 
extended action sum S among the variations in Vd and such that the solution 
sequence belongs to Cd- 

Otherwise said, we seek discrete paths c G Cd such that dS(c) ■ 5c = 0, for all 
5c G V d , with 6c(0) = = 5c(N). Using Theorem^ we get 

N-l 

= dS(c) ■ 5c = ^ D T>ELLd(tk-x,qk-i,tk,qk,tk+i,<lk+i) • (5t k ,5q k ) , 

k=l 

for all (8tk, 5q k ) G T>n k qk \, 1 < k < N — 1. Hence, the extended discrete Lagrange- 
d'Alembert (EDLA) equations read 

f DuEhLd(tk-i,qk-i,tk,qk,tk+i,qk+i) G ^ th ,q k ) ' 1 < k < N — 1 

\ (tk, Ik, tk+li qk+i) £ T^d, 

where T>° denotes the annihilator of T>. Let lo^ : Q x Q — > R, aG {1,... , m}, be 
smooth functions whose annihilation defines locally T>d, and let cu a : Q — > 
a G {1, . . . , m} be one- forms on Q locally defining 2?°. Since T*(T>) = TR, the latter 
ones are of the form u>(t,q) = (0,uj(t,q)), where with a slight abuse of notation we 
denote in the same way the component of the one- form in T*Q and the one- form 
itself. The EDLA equations can then be written as, 

DiLd(t k ,qk, tk+i,qk+i) + D 3 Ld(tk-i,qk-i,tk,qk) = , 
D 2 L d (t k ,q k ,t k+1 ,q k+ i) + D4,L d (tk-i,qk-i,tk,qk) = K^> a (tk,qk), (6) 
w d(^fc> tk+i, qk+i) = . 
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Notice that the discrete Lagrange-d'Alembert principle is not truly variational, in 
the sense that it does not correspond to the extremization of any action sum. This 
is in accordance with the nature of its continuous counterpart. Alternatively, we 
will refer to the EDLA algorithm © as a nonholonomic integrator, by analogy with 
the unconstrained case. 

Remark 2 (Well-posedness of the discrete problem). As it is also the case in 
unconstrained discrete mechanics |13j . the existence of solutions for the extended 
equations is not always guaranteed. If the mapping 

DiXT^ T*Q 

(to, <fo> h,q%, A) h-> (t , q , D 1 L d (t , q , tx,qx), -D 2 L d (t ,q , t 1 ,q 1 ) + X a uj a (t , q )) , 

is a local diffeomorphism, then for a pair (tfe_i, qk-x)i (tk, qk), there exists (ifc+i, qk+i) 
verifying the EDLA equations 10. The problem now arises from the fact that 
tk + \ > tk is not guaranteed, and therefore one might obtain inconsistent solutions. 
Nevertheless, one can ensure that, for specific choices of discrete Lagrangians |13| 
of natural (kinetic minus potential energy) systems, this situation does not occur 
away from points where the discrete energy is near zero. 

In the remainder of the paper, we assume that the EDLA equations © are well- 
posed and therefore induce an extended discrete Lagrange-d'Alembert map $ : Q x 
Q — > Q x Q, (t k -i,qk-i,tk,qk) (*fe)?fc)*fc+i)9fc+i)- The actual implementation 
of the EDLA algorithm can be carried out building on the discussion in [51 1131 

4. Geometric properties. In this section, we examine the geometric properties 
of the integrators derived from the discrete Lagrange-d'Alembert principle proposed 
above. It is important to keep in mind that the continuous flow of a nonholonomic 
Lagrangian problem does not have the same properties as the unconstrained flow |7|: 
on the one hand, the Poincare-Cartan form CIl is no longer preserved in general. 
On the other hand, the action of a symmetry Lie group does not generally give rise 
to momentum conserved quantities. However, the nonholonomic flow does enjoy 
some nice geometric properties with respect to these objects, which will guide our 
study of the corresponding discrete mechanics. 

Symplectic form: Consider the restricted action S : Q x Q — > K, 

S(t ,q ,tx,qx) = S(c) , 

where c is the unique solution of the EDLA equations satisfying c(0) = (io)(?o)> 
c(l) = (ti,qi). Using Theorem |21 with N = 2, we compute 

dS = \ a co a (tx,qx)+$*ei d -Ql d , 

and therefore conclude that $*f2£, d = ^L d + dfid, with (3d G T>°. Note that this 
is the discrete version of the behavior of the nonautonomous continuous flow with 
respect to the Poincare-Cartan two-form, Cx^l = d(3, with j3 E (2? u )° (see 0). 

Nonholonomic momentum map: Assume that the extended discrete non- 
holonomic Lagrangian system {L$, T>) is invariant under the (diagonal) action 
of a Lie group G on Q, that is, all the three elements are G-invariant. Let V denote 
the bundle of vertical vectors with respect to the canonical projection ir : Q — > Q/G, 

%?) = I£ e 0} • 

Among these symmetry directions, we are interested in selecting those ones which 
are also compatible with the nonholonomic constraints, that is, 

9(t, q ) = e q \ £Q(t, q ) e V iM) } . 
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Let 0° denote the (generalized) vector bundle over Q whose fiber at (t, q) is given 
by 9(t,q)- We now define the discrete nonholonomic momentum map as, 

Jf: QxQ — » ( 2 V )* 

(t ,qo,h,qi) i— > J2 h (to,qo,h,qi) : Q° -► R 

£ i-> (Jd(*o,go,*i,gi),0 • 

Note that this mapping is just the restriction of the usual discrete momentum map 
to the fiber bundle g v . Now, take a C°°-section of g v — > Q, that is, a mapping 
£ which for each (t, g) G Q gives us a symmetry direction £(t, (?) whose associated 
fundamental vector field lies in the constraint distribution. 

Proposition 1. Assume that (Ld,T>d,'D) is invariant under the action ofG. Then, 
the discrete time evolution of the nonholonomic momentum map is governed by the 
discrete momentum equation, 

(Jf(t 1 ,q 1 ,t 2 ,q 2 ),0 - (J? fc (to,gb,tl,?l),£) 

= (0+ (t U q 1: t 2 , q 2 ), (£{t 2 , q 2 ) - £(h, qi)) Q (t 2 , q 2 )) • (7) 

Proof. The Lie group G acts on Cd by means of the pointwise action. Then, 

N—l 

(dS(c),Zc d (c)) = E^'%xq) = 0- 

k=0 

On the other hand, since the space T>d is G-invariant, Cd is preserved by the group 
action. All this, together with the invariance of T>, implies that the solutions to the 
EDLA equations 10 are preserved by G, i.e., <& o ty g = \& g o $. 

Let (to, qo, ti, qi) G Q x Q and consider the corresponding solution to the EDLA 
equations. Take N = 2 and then, 

= (dS(c),£c d (c)) = (dS(t ,qo,h,qi),^Q X Q(t 0> qo,ti,qi)) 

= (x a u a {ti,qi) + $*e+ d -ei d ,^Q X Q(t , qo , h, qi )) . 

Now, if £g(ti,<7i) belongs to P(t ligi ), we deduce that 

( e L <i ( f l'9l' t 2,g2),CQ(^2,92)} = (0Z d (*O,?O,*l,«l),$Q(*O,?o)> ■ 

Finally, 

{J2 h (h, qi ,t 2 ,q 2 ), - (J2 h (to,Qo,h, qi ),0 

= ( e i d (*l)5l)*2,Q52), (£,{t 2 ,q 2 ))-Q(t 2 ,q 2 )) - (@l d (t ,q ,ti,qi), (£(ti, ?i))q{*o, go)) 

= ( Q i d {h,gi,h,Q2), (£(t 2 ,q 2 ))-Q(t 2 ,q 2 ))-(®L d {tx,qi ) t 2) q 2 ), 9i))q(*2, 92)) , 

which is the desired result. □ 

A distinguished class of sections of the bundle fl 77 is formed by the constant 
ones, £(g) = £. They correspond to elements £ of the Lie algebra which always are 
compatible with the constraints, that is, £g(t, q) G 2?(t. q ), for all (t,q) G Q. These 
special elements are called horizontal symmetries in the literature of nonholonomic 
mechanics 0O3E]. 

Corollary 1. If £ G Q is a horizontal symmetry, then the associated component of 
the discrete nonholonomic momentum is preserved by the EDLA algorithm. 

Proof. It is immediate from {7}, since in this case £(t 2 , q 2 ) — 51) = 0, and hence 
(JS h (h,q u t 2 ,q 2 ),0 = (JS h (to,q ,ti,qi),0. □ 
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Nonholonomic Chaplygin systems: It may also happen that the generalized 
bundle g v over Q is trivial, that is, 0u,q) — f° r au (*> 0) G Q- I n this case, there is 
no nonholonomic momentum map and hence we must look for different geometric 
properties of the flow other than Proposition ^ Under the additional hypothesis 
D + V = TQ (dimensional assumption, cf. PIE]), we deduce that T> complements 

V in the tangent bundle of Q, and therefore constitutes the horizontal space of a 
principal connection. We denote its associated connection one-form by A : TQ — » g. 
This class of nonholonomic systems are called generalized Chaplygin systems . It 
is known that, after the reduction by the action of the Lie group, these systems give 
rise to an unconstrained system subject to an external force of gyroscopic type. In 
the following, we show that the discrete mechanics also shares this feature. 

Assume that the discrete constraint space T>d and the action are such that TT> c i + 

V x V — TQ (an hypothesis that we term discrete dimensional assumption). Let 
7r : Q — > Q/G be the canonical projection, and consider the map 

v : V d /G — > Q/G x Q/G 

[(to,qo,h,qi)] i — > (n(to, qo),w(ti, qi)) ■ 

Note that both spaces have the same dimension due to the definition of T>d and 
the dimensional assumption. Indeed, diml?^ = dim I? = 2dimQ — dimG. On the 
other hand, if p : Vd — > Vd/G denotes the projection from T>d to its reduced space, 
then one can verify that kerp* C TT>d n (V x V). By a dimensional argument, we 
conclude that kerp* = TT>d fl (V x V), and therefore v is a local diffeomorphism. 

We say that T>d is right-rigid with respect to the G-action ip if the following 
property holds: given (t ,q ,ti,qt) G V d and g G G, if (t ,q ,g(ti,qi)) G V d , 
then g = e (where we are using the abbreviated notation g{t\,qi) — ip(g,ti,qi))). 
Clearly, if T>d is right-rigid and invariant under the diagonal action, it is also left- 
rigid. Intuitively, the notion of right-rigidity (resp. left-rigidity) means that T>d is 
not invariant under the action Id xtp : GxQxQ —> QxQ, (Id xijj)(g, to, qo, t\, q±) = 
(t ,q ,ip(g,ti,qi)) (resp. V x Id). 

Proposition 2. Let (Ld,T>d,T>) be G- invariant. Assume that the generalized bundle 
q v on Q is trivial and that the discrete dimensional assumption holds. Then, if T>d 
is right-rigid, the local diffeomorphism v is global. 

Proof. We use the abbreviated notation q = (t,q) G Q. Take [(5oi9i)L [(<7oi<7i)] e 
T>d/G such that v[(q ,qi)] = f[(<7o>?i)]' Then there exist go, gi G G such that q Q = 
3o9o> Qi = 9ili- Since (3o><Zi) G V d , then (go9 ;5o9i) G V d by G-invariance. Alter- 
natively, we have (g , q x ) G T> d and, at the same time, (go%^ 9oQi) — (<Zoj 9o9i l qi) G 
V d . Now, by rigidity, we conclude g = g x , and hence [(g >9i)] = [(<?o>ti)]- □ 

Therefore, under the global identification provided by v, we can define a reduced 
discrete Lagrangian L* d : Q/G x Q/G -> K, L* d (Tk,rk+i) = ^(^ _1 (f fc , r fc+ i)), where 

: (Q x Q)/G — > R is the reduction of L<j to (Q xQ)/G, and we regard T>d/G as a 
submanifold of (Q x Q)/G. Locally, if we identify Q with Q/G x G, (t, q) — (r,g), 
then we can take local coordinates (ro,Fi, /o,i) G Q/GxQ/GxG on (QxQ)/G. In 
this way, the projection Q x Q —> (Q x Q)/G reads (ro,go,^i,3i) l— * (^Oj^1j/o,i = 
g ~ 1 gi). Moreover, when regarding T>d/G as contained in (Q x Q)/G, we have that 
/o,i = /o,i(^07^i) for (^Oj^ii/o.i) G T>d/G. Finally, if the G-action acts trivially on 
the time component of Q, we can further write Q/G = K x Q/G, r = (t,r). Now, 
we are in a position to state the following result. 
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Proposition 3. Under the hypothesis of Prop. \^ assume ip acts trivially on the 
time component of Q. Then, the solutions of the ED LA equations project onto the 
solutions of the reduced extended discrete Lagrange-d'Alembert (REDLA) equations, 

DiLj(tk,r k ,tk+i,r k +i) + D 3 Lj(t k -i,r k -i,tk,r k ) = 0, 

D 2 L* d (t k , r k , t k+ i, r k+ i) + D 4 L* d (t fc _ i , r fc _ i , t k , r k ) = F~ (t k ,r k ,t k+1 , r k +i) 

+F + (t k ^ 1 ,r k ^ 1 ,t k ,r k ) , 

where the expression of the forces in bundle coordinates is given by 

F-(r k ,r k+1 ) = g r fc fc +1 ( T k,r k+1 ) - Rf k , k+1 (r k ,r k+1 ){Ai oc (r k )(-))^ , 

F + (r k -i,T k ) = Q ^ - ( d ^' k (rk-i,f k ) + L fk _ lik( r h _ u r h ){A loc {r k )(-)) \ , 

where Ai oc (r) is the local form of the connection one-form A. 

This result can be proved using a similar argument to the one carried out in [5] for 
autonomous systems. The REDLA integrator is an appropriate version for nonau- 
tonomous systems of the generalized variational integrators developed for systems 
subject to external forcing in |14j . This is in accordance with the situation in the 
continuous case where, as we mentioned before, the reduction of the Chaplygin 
system gives rise to an unconstrained system subject to a gyroscopic external force. 

Energy conservation for autonomous constrained Lagrangian systems: 
The discrete system (L dl 'D d ,D) is autonomous if L d , T> and T> d are invariant under 
the additive action of R on the time component of Q. In this case, q v = g = R. 
Consequently, the discrete nonholonomic momentum map coincides with J d , which, 
as we have already seen, is given by JL d (to,Qo,ti,Qi) = — (to, 1o, ti, qi)dti = 
—E^ d {to,qo,ti,qi)dt^. Corollary ^ thus yields 

E^ d (t k ,q k ,t k+1 ,q k+ i) = E Ld (tk-i,qk-l,tk,Qk) ■ 

Proposition 4. If the discrete system (L d ,T> d ,T>) is autonomous, then the EDLA 
algorithm preserves its associated discrete energy. 

5. Conclusions. We have proposed a discrete version of the Lagrange-d'Alembert 
principle for nonautonomous Lagrangian systems with nonholonomic constraints. 
We have studied the geometric properties of the integrators derived from this prin- 
ciple, paying special attention to the evolution of the symplectic form and the 
nonholonomic momentum map, and the conservation of energy. Future work will 
be devoted to develop a numerical error analysis of these integrators making use of 
backward error techniques. 
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